Spin-Charge separation in a model of two coupled chains 
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Q^ . A model of interacting electrons living on two chains coupled by a transverse 

hopping t_L, is solved exactly by bosonization technique. It is shown that t±_ does 
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\^ . modify the shape of the Fermi surface also in presence of interaction, although charge 
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short distances, x <C ^ = {up — Mo-)/4t_L, the two chain model is not sensitive to 
t_L, while for larger separation x ^ ^ inter-chain hopping is relevant and generates 
further singularities in the electron Green function besides those due to spin-charge 
decoupling. 
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A central problem in the theory of strongly correlated systems is whether the mechanism 
leading to the breakdown of Fermi liquid (FL) theory in one dimensional (ID) models can 
be generalized to higher dimensions [jl|. In most ID systems, the non FL behavior manifests 
itself in two distinct ways: the correlation functions show power law behavior with coupling 
dependent anomalous exponents [0] leading, in particular, to the smoothing of the Fermi 
surface that is often considered as the key signature of non FL behavior. An independent, 
but equally effective mechanism which induces the breakdown of FL without generating 
anomalous exponents, is the well known spin-charge decoupling [^ which originates from 
the dynamical independence of charge and spin excitations. Many unsuccessful attempts 
have been made in the past years to search for a breakdown of FL theory in 2D through 
the study of the discontinuity in the momentum distribution, i.e. by looking for anomalous 
exponents. In this Letter we address the problem of the breakdown of FL in 2D models 
0J^ by looking directly for spin-charge decoupling. In particular, we analyse the stability 
of spin-charge decoupling with respect to the introduction of transverse hopping between 
two chains. We formulate and solve via bosonization technique, a model which shows non 
FL behavior without anomalous correlation exponents and we calculate several physical 
quantities. At low energy we find new collective excitations related to the transfer of an 
electron across the chains which modify the analytic structure of the Green function without 
restoring FL behavior. Some results of the two chain problem can be easily extended to an 
array of A^-chains, i.e to strongly anisotropic 2D systems. 

A simple ID model which shows spin-charge decoupling without anomalous power-law 
decay in the correlation functions, can be obtained by including only forward scattering 
processes between electrons on the same branch of the Fermi surface (usually referred to as 
gi interactions). The hamiltonian of this model is [Q: 

Ho = Yl vpik - kp)al^a^^ + ^ vpi-k - kF)bl^b^^ + 

fc,(T k,a 
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where the operator a[^(&[^) creates an electron of momentum k and spin a belonging to the 
branch with positive(negative) slope vp{—vf)- Due to the linear dispersion relation of the 
hamiltonian ^ (|T]), this model can be solved by use of standard bosonization technique being 
a special case of the well known Tomonaga-Luttinger model. While this interaction does not 
affect the ground state properties, the excitations are profoundly different. Two low-lying 
gapless excitations with linear dispersion, are found, describing charge and spin collective 
modes with different sound velocities Up = Vf + {g^\\ +5'4±)/2vr and u^ = Vf + {g4\\ — 5'4±)/27r. 
No quasi-particle excitation exists and as a consequence the Green function does not have 
simple poles but branch cut singularities. The bosonization method allows to calculate the 
exact form of the Green function in real space and time: 
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{x + Upt — ir/signt) {x + u^t — irjsignt) 
After Fourier transforming, the spectral weight, e.g. for the right moving electrons, can be 
easily obtained (for more general ID models see Ref. 7): 

[e{q)e{uj - u^q)e{Upq - uj) + e{-q)e{uj - Upq)e{u^q - uj)] 

A[Kp + q,uj) — , • [o) 

n J\u — Ufjq\\uj — Upq\ 

In the particular case g^w = g4± = the spin and charge velocities u^ and Up coincide 

and the branch cut merges in a simple pole reproducing the standard free particle Green 

function. In this simple model, although the low lying excitations can not be described in 

terms of the Landau Fermi liquid theory, all the equal time correlation functions at zero 

temperature, including the momentum distribution, coincide with the non-interacting ones. 

The form of the spectral function (^) explicitly shows that non FL behavior can occur even 



without anomalous exponents. In this case, the breakdown of FL is related to a change in 
the analytic structure of A{k,uj) which does not generate the logarithmic singularities in 
perturbation theory present in many ID models. 

This model is a prototype to study the relevance of spin-charge decoupling in a system 
of chains coupled by transverse hopping. As a first step let us consider the simple case of 
two chains. The total hamiltonian will be the sum of two terms like Eq. (|Tp, one for each 
chain, plus an hopping term between the two chains: 

H± = -t± J2 («L,i«fca,2 + bL,Aa,2 + H.c.) (4) 

kcr 

where the suffixes 1 or 2 refer to the two chains. We will now show how it is possible to solve 
exactly this problem by bosonization technique. As there is no term in the Hamiltonian 
which couples right to left moving electrons, we treat only the case of the right moving 
electrons. The hamiltonian Eq. ([l|) can be easily rewritten in terms of the density operators 
of the electrons on the positive branch : 

27r 

Ho = ^Mpl] [pi(g)pi(-g) +P2(g)p2(-g)] + 
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+ -r'^'^ ^ Wi{q)M-<i) + ^2(g)o^2(-g)] (5) 

^ q>0 

where Pi(?)(c"i(g)) is the charge(spin) density operator of the right moving electrons on chain 
1. In order to write the transverse hopping Eq. (^ in terms of the density operators, it is 
necessary to introduce the boson representation of fermion operators [Q . Then one obtains: 

H± = 2it_L- — / dx [sin {(f)ps{x) + (pasix)) + sin (0p^(x) - (pasix))] (6) 

where the operator 0j's are the phase fields related to the densities: 

(7) 
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and they are defined as [§,|1: 









The suffix i labels any of the four densities in Eq. (J^), and the factor exp(— ag/2) plays 
the role of an ultraviolet cutoff for divergent integrals. In terms of the densities (|^), the 
hamiltonian (^) remains diagonal. Let us introduce the fermionic fields corresponding to 
the densities (0). Then, the hopping term (|) acquires a simple form: 

-t^Jdx {^l{x)^^,ix) + H.c.) + {e'''-^<iJl{x)<ijUx) + H.c.) . (8) 

By defining the Fourier transform Cp{k) {ca{k)) of the field operator \l/ps {"^as), the total 
hamiltonian, i.e. Eqs. (|^) plus (H), can be rewritten as: 

H = T. {Mk - kF)cl{k)Cp{k) + u„{k - kF)cl{k)cM) + 

k 

-^±E(ci(^)c.(^) + H.c.) + 

k 

-t±Yl (cp(^ + '^F)ci{k - kp) + H.C. + cl{k - kF)cl{k + kp) + H.c.) . (9) 

fc>0 

This hamiltonian is the sum of bilinear terms of fermion operators, therefore it can be easily 
diagonalized. By performing a particle-hole transformation and then a unitary transforma- 
tion, four excitation branches are obtained: 

ei(?) =Upq 
62(9) = u„q 
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t^{q) = \{up + u^)q + y[l{up - u^)q) + W]_ 
e4(g) = \{up + u^)q - \j{^\{up - u„)q^ + U\ 
where q = k—kp is positive. Having set the chemical potential equal to zero, the ground state 

is obtained by filhng the branch 64 (g) up to a momentum Q = 2t±/ ^UpUa^ corresponding to 

e4(Q) = 0. 

We have checked that, although this spinless fermion excitation spectrum seems quite 
complicated, it correctly reproduces all known results, in particular the charge- and spin- 
density normal modes in the trivial limits a) g^w = = g4^± for t± 7^ 0; b) g^w 7^ 0, g4± 7^ 
for 1 1 = 0. 



From the ground state and the excitation spectrum several physical properties of the 
system can be calculated. The ground state energy is given by twice the integral of 64 (g) 
from g = up to g = Q- This accounts for both right and left moving electrons, giving: 

AE 1 4ti u 

L 27r [Up - Ua) Ua 

Notice that Eq. (|TT|) reproduces the correct result in the non-interacting case, i.e. when 
Up = Ufj = vp. The energy correction ([TI|) being proportional to t]_ implies a non zero value 
of the transverse hopping operator averaged on the ground state, which in turns means 
that spin-charge decoupling is not sufficient to generate confinement of the electrons within 
each chain |Q. In order to analyse this issue more deeply, we calculate the difference of 
the occupation numbers between the bonding and the anti-bonding band. By labelling the 
two bands according to the correspondent transverse momenta, i.e. k± = Q for the bonding 
combination, and k_j_ = n for the anti-bonding one, we get: 

L 2iT Up- u^ \u^J 

which, again, does not show confinement in the two chain problem. As a further probe 
for confinement, we have evaluated the number fiuctuations between the two chains. This 
quantity is easily related to the long wavelength limit of the density structure factor, which 
can be calculated by the bosonization method: 

2 L n{up-u^) \ ^ Up) ' 

This expression, apart from the trivial dependence upon t±, is a function of the interaction 
through Up and u^. In the simple case of spin-isotropic interaction, i.e. g4\\ = g4±, the 
number fiuctuation is a monotonic decreasing function of the interaction, which goes to zero 
only when the coupling tends to infinity. 

It is now interesting to investigate the stabihty of spin-charge decoupling with respect to 
the introduction of transverse hopping. From a simple inspection of the excitation spectrum 
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in Eq. (0), one can argue that for large momentum the simple ID picture is recovered, 
together with spin-charge decoupling. In fact, for large q, e^i^q) — > ei(g) = Upq, e^i^q) — >■ 
€2(q) = u^q, and the spectrum tends to the t_L = case. On the other hand, new excitations 
appear at low energy, centered around q = 2t±/y/UpUa, with an almost linear spectrum and 
velocity 2upUcr/{up + Ua)- The length scale which separates the two regimes is given by 

^ = {up-u^)/4t^. (14) 

This form of the excitation spectrum leads to important consequences in the spin and 
charge density-density correlation functions. On chain 1 the latter is defined as 

-i^(t)([pi(g,t),pi(-g)]) = -^{([p,(g,t),p,(-g)]) + ([p,(g,t),p,(-g)])} 

The first term on the right hand side is unaffected by t_j_ and it contributes to the spectral 
weight A{uj,q) with a delta-function centered at u = —Upq. For q^ ^ 1, the t_L ~ 
regime is recovered and most of the spectral weight lies in the sharp peak at cj ~ —Upq 
as in the case of two independent chains (Fig. 1). In the opposite limit g^ <^ 1 the effect 
of transverse hopping dominates over interaction and the system behaves like in the non 
interacting case where the ps term contributes to the spectral weight with two delta-peaks 
centered respectively at cj = —vpq + 2t_L and aX uj = —vpq — 2t_L (see Fig. 2). In both 
limits, the interplay between the transverse hopping and the interaction, has the effect of 
broadening the peaks in A{k,uj), which now acquire a finite width. 

A direct probe of spin-charge decoupling can be obtained through the single-particle 
Green function. This property is difficult to extract from our solution, because the fermionic 
operators have a very complicated expression in terms of the normal modes which diago- 
nalize the hamiltonian. However, it is possible to obtain its asymptotic behavior at large 
distance which contains the relevant information regarding FL behavior. In this regime, we 
approximate the excitation branch €4(5) in Eq. ( PUJ ) near q = Q with a linear spectrum and 



velocity Ur = 2upUa/{up + u^)- This simplification allows to calculate the real space Green 
function at fixed transverse momentum {k± = 0, vr): 

{^oix,t)^l{0, 0)) ~ ei(kp+^kp)-^x - Upt)-I{x - u^t)~l{x - M.t)-^ (15) 

where 

Akp = — log (up/u^). (16) 

Up - Ua 

An analogous expression can be derived for the anti-bonding combination of the two chains 
{k± = it) by reversing the sign of Akp. Equation ( [T6| ) shows that the transverse hopping 
does move the Fermi points of the bonding with respect to the anti-bonding band. The 
asymptotic expression Eq. ( p^f ) corresponds to the region x, vpt ^ C,', for smaller separation 
we expect to recover the single chain limit Eq. @. Therefore, we find a crossover between 
two regimes. On a scale x < ^ fluctuations between the two chains are suppressed and 
the model essentially behaves like two independent chains. At long distance (x > .^), the 
coupling between the chains is reintroduced but interaction still plays a crucial role. An 
electron excitation now decays into a triplet of elementary excitations, and not just a holon 
and a spinon. 

Encouraged by the absence of singularities in the transverse hopping as it comes out 
from the exact solution, we checked the bosonization results against perturbation theory in 
t±. The first order correction to the single particle Green function can be easily obtained 
up to first order in t±: 

Go w ik,u) = G^o) ik,u)Tti. [g(°) ik,cu)f (17) 

where the — (+) corresponds to k± = 0{n), and the unperturbed Green function G^'^\k,uj) 
is just the Fourier transform of Eq. (§). From Eq. (|l^) we can calculate the correction to 
the bare momentum distribution Uq (k) = n^\k) = 9(kp — k): 
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This integral gives apparently zero due to a well known anomaly of the T = perturbation 
theory [|l0|. This problem can be avoided by working at finite temperature. In this way the 



momentum distribution becomes: 

(5no(/cF + q) = t±-, ^ [fi(3uaq) - fiPupq)] 

[Up - Ua)q 

where the function f{x) is the Fermi distribution. By taking the T — > limit, we get: 

6no{k) = log {up/ua)6{k - kp) 

Up -u„ 

which can be interpreted as a shift of Fermi momenta given by Eq. (|16|). Analogously, we 
can easily check that the total energy correction induced by this shift coincides with Eq. 
(p!T|). Therefore, at least up to lowest order in t±, perturbation theory agrees with the exact 
solution. 

It is straightforward to generalize these perturbative results to the case of an array 
of A^ chains. The relevant equations remain unchanged, the only difference being that 
the inter-chain hopping operator acting at transverse momentum k±_ is formally replaced 
by —t±cos{ki_). For example the correction to the longitudinal Fermi momentum with 
transverse momentum k±_ is: 

A , /, N t_LCOs(/c_L) , / / N 

Afc^(fcx) = ^^log(M^K)- 

Up -Ua 

to first order in t^. This equation gives the shape of the Fermi surface for a strongly 
anisotropic 2D system as a function of the spin and charge velocities. However, we have 
not been able to find the exact solution of the A^-chain problem to all orders in t±_ and the 
crucial issue of the breakdown of FL behavior in two dimensions is still open. 
We thank E. Tosatti, J. Voit, and P. Nozieres for helpful discussions. 
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FIGURES 
FIG. 1. Charge density-density spectral function in the channel ps = l/v2(pi — P2)- The 
calculation has been performed setting Up — u^ = 0.1, t±_ = 0.01 and q = 0.7 in units of vp = 1. 

FIG. 2. The same as Fig. 1 but at momentum q = 0.05. 
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